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Abstract 

Due to its intimate relation to Spectral Theory and Schrodinger op- 
erators, the multivariate moment problem has been a subject of many 
researches, so far without essential success (if one tries to compare with 
the one-dimensional case). In the present paper we reconsider a basic 
axiom of the standard approach - the positivity of the measure. We in- 
troduce the so-called pseudopositive measures instead. One of our main 
achievements is the solution of the moment problem in the class of the 
pseudopositive measures. A measure fi is called pseudopositive if its 
Laplace-Fourier coefficients [ik,l (r) , r > 0, in the expansion in spherical 
harmonics are non-negative. Another main profit of our approach is that 
for pseudopositive measures we may develop efficient "cubature formu- 
las" by generalizing the classical procedure of Gauss-Jacobi: for every 
integer p > 1 we construct a new pseudopositive measure v v having 
"minimal support" and such that fi(h) = v v (h) for every polynomial h 
with A 2p h = 0. The proof of this result requires application of the famous 
theory of Chebyshev, Markov, Stieltjes, Krein for extremal properties of 
the Gauss-Jacobi measure, by employing the classical orthogonal polyno- 
mials pk,i;j, j > 0, with respect to every measure fik t i- As a byproduct we 
obtain a notion of multivariate orthogonality defined by the polynomials 
Pk,i-,j- A major motivation for our investigation has been the further devel- 
opment of new models for the multivariate Schrodinger operators, which 
generalize the classical result of M. Stone saying that the one-dimensional 
orthogonal polynomials represent a model for the self-adjoint operators 
with simple spectrum. 

1 Introduction 

The univariate moment problem is one of the cornerstones of Mathematical 
Analysis where several areas of Pure and Applied Mathematics meet - contin- 
ued fractions, quadrature formulas, orthogonal polynomials, analytic functions, 
finite differences, operator and spectral theory, scattering theory and inverse 
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problems, probability theory, and last but not least, control theory, see e.g. the 
collection of surveys in [23 and the comprehensive recent account 03] on the 
numerous applications of the moment problem to spectral theory. On the other 
hand, the multivariate case is much more complicated, and we refer to fT5). 
[37], 0U|, 0U and the references given there for some recent developments. 
However, the state of the art in the multivariate moment problem seems to be 
well characterized by a remark in the versatile survey |20l p. 47], saying that 
only comparatively little from the comprehensive theory of the classical moment 
problem has been extended to dimension d > 1. 

The main purpose of the present paper is to introduce a modified moment 
problem for which the solutions are in general signed measures and belong to 
the class of what we call pseudo-positive measures. The motivation for this new 
notion is the possibility to generalize the univariate GauB-Jacobi quadratures 
to the multivariate setting, and thus to approximate multivariate integrals in a 
new stable way. Let us emphasize that we do not claim to solve the multivariate 
moment problem in its classical formulation. 

In order to make our approach clear, let us first recall the usual formulation 
of the multivariate moment problem: it asks for conditions on a sequence of real 
numbers c = {c Q } QgN d (here N denotes the set of all non-negative integers and 
we use the multi-index notation x a — x^x^ 2 • ■ ■ %d d ), such that there exists a 
non-negative measure \i on M. d with 

c a = / x a dfi (x) (1) 

for all a £ Nq. Let us denote by C \x\, X2, Xd] the space of all polynomials 
in d variables with complex coefficients. With the sequence c = {c a } QeN <i we 
associate by linear extension to C [x\, x^, Xd] a functional T c , by putting 

T c (x a ):=c a for q e N d . (2) 

By a theorem of Haviland, a necessary and sufficient condition for the existence 
of a non-negative measure \x satisfying Q is the positivity of the sequence c — 
{c Q } QGN g, i.e. P > implies T c (P) > for all P £ C [x 1 ,x 2 , ...,x d ] (here P > 

means that P (x) > for all x £ R d ), cf. 7, p. 111]. As is well known, 
the condition of positivity is difficult to apply, and in practice one uses the 
weaker and easier to check algebraic condition that the sequence c = {c Q } a(EN d 
is positive definite; by definition the sequence {c a } a ^d is positive definite 1 if 
and only if 

T C (P*P) > for all P £ C [xi,x a , -,x d ] ; 

here P* is the polynomial whose coefficients are the complex conjugates of 
the coefficients of P. Let us remind that for d = 1, positivity and positive- 
definiteness of T c are equivalent. However, for d > 1 this is not true, and this 

1 Some authors call such sequences " positive" , while we are closer to the terminology of 
say 0. 
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is a consequence of the fact, already known to D. Hilbert, that there exist non- 
negative polynomials which are not sums of squares of other polynomials, see 
e.g. [221, 0. 

As mentioned above, we shall consider a different setting of the moment prob- 
lem. Let us postpone at the moment the motivation for our approach, and let us 
concentrate on our new setting which needs some technical preparations from the 
theory of harmonic functions. We assume that for each k = 0, 1, 2, the func- 
tions YJy : M. d — ► R, I = 1, ak, form a basis of the set of all harmonic homo- 
geneous complex- valued polynomials 2 of degree k € No, and they are orthonor- 
mal with respect to the scalar product (f,g} sd _ 1 := L d -i f (@) 9 {6)d9, where 
^ d ~ 1 := {x e M d : |x| = 1} is the unit sphere, and r = \x\ = ^Jx\ + .... + x 2 d 
is the euclidean norm, cf. [H] or For x £ M. d we will use further the 

representation 

x = r6 for r > 0, 6 e § rf_1 . 

The functions Y^j are called solid harmonics and their restrictions to 
spherical harmonics. An important property of the system | | 2j Y^j (x) , j,k G 
No, I = l,...,afc, is that it forms a basis for C [xi, X2, Xd\- This result fol- 
lows from the GauB decomposition of polynomials which says (cf. [SJ Theorem 
5.6, Theorem 5.21, p. 77 and p. 90], gt], or |H1 Theorem 10.2]) that every 
polynomial P may be expanded in the following way, 

[dcg(P)/2] + l 

P(x)= ]T \x\ 2j Pj (x), (3) 

3=0 

where pj are harmonic polynomials, degP denotes the degree of P and [x] 
is the integer part of a real number x. Since each pj is a linear combination 
of the solid harmonics Y^i (x) , k € No, I = 1,2, ...,<!&, it is clear that the 
system | a; | 2j Ykj, (x) , j, k € No, I = 1, ctfe, is a basis for C [x%, X2, Xd\- It is 
instructive to discuss the relationship between the GauB decomposition and the 
Laplace-Fourier series: recall that for a sufficiently nice function / : W 1 — * C 
(e.g. continuous) the expansion 

oo a k 

/w=EEa,iW^(»), (4) 

fe=0 1=1 

is the Laplace-Fourier series with the Laplace-Fourier coefficients given by 

h,i{r)= [ f{rO)Y k j{6)d6. (5) 

Suppose now that / is a polynomial: then © implies that the Laplace-Fourier 
series 10} is a finite series, and the functions fk,i (r) r~ k are polynomials in the 

2 For a reader not familiar with the spherical harmonics it will be enough to consider the 
two-dimensional case d = 2 where the basis is simple and given in Section 18.31 On the other 
hand in Section |7| we develop our theory in the case of the strip where all what one needs is 
expansion in Fourier series. 
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variable r 2 . Moreover from @ directly follows that each fk,i (r)r k , k 6 No, I = 
1, 2, cifc, is a polynomial of degree < 2s — 2 if and only if for all x € R n 

A s /(.t) = 0. (6) 

Here A denotes the Laplace operator defined by A = + ... + -g^, and A s is 
the s-th iterate of A for integers s > 1 . In view of © let us recall that a function 
/ defined on an open subset U in R d is polyharmonic of order s if A s / (x) = 
for all x eU, see [3]. 3 

Now we come to the cornerstone of our approach. Let T c : C [xi, x<i, a; J — > 
C be a functional associated to a sequence of moments c a ,a S Nq. Using the 
basis | | 2-3 Ykj (x) : j,k £ No, I — 1, at one can define a problem equivalent 
to the usual one Q, by means of the sequence {Cj"}j£$ defined as follows, 

cf' l) := T c (|x| 2j F M (x)) for j = 0, 1, 2, ... (7) 

In order to distinguish them from the usual moments c a , the numbers Cj are 
sometimes called distributed moments, see ^0], |H| - |25j . |2fi| . |27j . |29j . 

We say that the sequence {c a } Q , gN <i or the associated functional T c is pseudo- 
positive definite if for every fixed pair of indices (k,l) with k G No and ^ = 

1, ...,<Zfc the sequences {c^ ^}jgN an d { c j+i}j'ePJo are positive definite. Equiv- 
alently, for each solid harmonic Y k i{x) ,k € No, I = l,...,afe, the component 
functional Tf. t i : C [xi] — > C defined by 

T k ,i (p) := T c (p{\x\ 2 )Y k ,i (x)) for p e C [x x ] (8) 

has the property that T*,; (p* {t)p{t)) > and T fe) ; (ip* (i)p(i)) > for all 

p e c [si] . 

Now we can formulate and successfully solve the following modified moment 
problem: Given a pseudo-positive definite sequence c = {c a } QeN d and its asso- 
ciated functional T c , find the conditions for the existence of a signed measure /j, 
on R d such that 

J P(x)dfx = T c (P) ior aH P £C[xx,x 2 ,...,x d ]. (9) 

Two remarks are important: first, we allow /i to be a signed measure on K", and 
this requirement is motivated by our constructive formulas for approximating 
integrals developed in later sections. Secondly, it follows from (j5J) that the 
measure /i, considered as a functional on C [xi, x%, xj[ , is pseudo-positive 
definite. The remarkable thing which will be seen from our further development 
is that problem JJjJ has a solution /i which is pseudo-positive which means that 
the inequality 

f h(\x\)Y kt i(x)dfj,(x)>Q (10) 

3 In the last section we provide some remarks on the significance of the polyharmonic 
functions in approximation theory which have motivated also the present research. 
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holds for every non-negative continuous function h : [0, oo) — * [0, oo) with com- 
pact support and for all pairs of indices (k, I) with k e No and 1 = 1,2, a k . 

As a first evidence that the pseudo-positivity is a reasonable generalization 
of the univariate positivity notion, we present in Section[21the following solution 
to the modified moment problem JjJJl, provided in two steps: 1. By a classical 
one-dimensional argument, for the component functionals T k j associated with 
the pseudo-positive definite functional T c : C [xi, xi, Xd] — * C there exist 
non-negative univariate representing measures fik,i on [0, oo) . 2. If they satisfy 
the summability assumption 

oo Qfc /.oo 

EE/ r N r- k dn ki i (r) < oo for all JVeN (11) 
k=o i=i •'° 

then there exists a pseudo-positive signed measure \i on M. d representing T c , i.e. 
J2J holds. Further the following important identity 

„ oo afc /.oo 

/ f(x)d^ = J2J2 fk,l(r)r- k dn k ,i(r) (12) 

JR n fe=Q l=l JO 

holds for any continuous, polynomially bounded function / : R" — > C; here 
fk,l i r ) are the Laplace-Fourier coefficients. Equation (jT2")l will be the key for 
defining polyharmonic Gaufi-Jacobi cubatures as we shall show below. 

Another strong supporting evidence for the nice properties of the notion 
of pseudo-positivity is the satisfactory solution of the question of determinacy 
in Section |3| we show that the representing measure fi of a pseudo-positive 
definite functional T : C [xi, X2, xj\ — > C is unique in the class of all pseudo- 
positive signed measures whenever each component functional Tkj defined in 
(jHJ has a unique representing measure on [0, oo) in the sense of Stieltjes (for the 
precise definition see Section And vice versa, if a pseudo-positive functional 
T is determinate in the class of all pseudo-positive signed measures and the 
summability condition l|ll|) is satisfied, then each functional Tkj is determinate 
in the sense of Stieltjes. The proof is essentially based on the properties of the 
Nevanlinna extremal measures. 

Let us illustrate the notion of pseudo-positivity in the case where the signed 
measure fi has a continuous density w (x) with respect to the Lebesgue measure 
dx. We put d/j, (x) = w (x) dx in 10 and take into account that Y^ ; i (x) — 
\x\ k Yk,i (0), and dx = r d ~ 1 d6dr, (for detailed computations see Proposition l41l) . 
We see that the component functionals T ki i defined in JSJ) are now given by 

poo 

T Kl (p) = / p (r 2 ) r k+d - l w k ,i (r) dr, (13) 
Jo 

where w kl i (r) are the Laplace-Fourier coefficients of the function w as defined 
in From (|13fl it is obvious that the non-negativity of w k .i implies that the 
measure \i is pseudo-positive and the corresponding functional T defined by © 
is pseudo-positive definite. We regard now d/j, k ,i (r) — r k+d ~ 1 w k j (r) dr as a 
univariate non- negative measure which represents the functional T k j. 
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Now we are moving to our main theme, the construction of GauB-Jacobi type 
cubatures for pseudo-positive measures. Let us first recall some terminology: 
By a cubature formula one usually means a linear functional of the form 



C (/) := aif (xi) + .... + a.f (x s ) = / (]T a,S (x - Xj ))f (x) dx (14) 



defined on the set C(]R Tl ), the set of all continuous complex- valued functions 
on R n ; here 5 is the Dirac delta function. The points xi, ...,x s are called nodes 
and the coefficients oil, a s £ R weights. A cubature formula C (•) is exact on 
a subspace U of C (W l ) with respect to a measure v if 

C(/) = J f{x)dv (15) 

holds for all / G U. If U s is the set of all polynomials of degree < s, and the 
cubature is exact on U s but not on U s +\, we say that C has order s. It is 
common to call a cubature in the case d = 1 a quadrature. 

In our construction we will use the Gaufi-Jacobi quadrature, so let us recall 
its definition: Let v be a non-negative measure on the interval [0, R] and s > 1 
be an integer. If the cardinality of the support of v is > s then there exist s 
different points tj in the interval (0, R) and s positive weights otj which define 
the classical Gaufi-Jacobi measure 

s 

dv^ {t) = ^ aj 5{t-tj)dt (16) 

3=1 

and the corresponding Gaufi-Jacobi quadrature C (/) = J ^ / (t) dz/ s ) (i) = 
2j=i a j/ fe) satisfies (|15H for all polynomials of degree < 2s — 1. For formal 
reasons we put if the cardinality of v is < s. 

It is not our intention to survey the numerous approaches to cubature formu- 
las; one may consult the references in |2H], JHj, 07], ^7J and in particular the 
monograph of S. L. Sobolev |44 | where the minimization of the error functional 
of the formula in (|15f) is the main objective. 

Our approach, which is rather different from the usual cubature formulas, 
is based on the identity Q12JI. In order to be precise, let us begin with the 
assumptions: let /i be a pseudo-positive measure and define the non-negative 
component measures [ik,i by the identity 

h{\x\)dii k ,i =j h(\x\)Y k>l (x)d» (17) 

valid for all continuous functions h : [0, oo) — > C with compact support. Let 
ip : [0, oo) [0, oo) be the transformation tp(t) = t 2 and let be the image 

4 The points tj are the zeros of the polynomial Q s (t) which is the s— th orthogonal with 
respect to the measure v on [0, R] . It is important that tj £ (0, R) if the support of v has 
cardinality > s, cf. Chapter 1, Theorem 5.2 in 1121 . and Theorem 5.1 in Chapter 3.5 in 1321 . 
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measure of (ik,l under tp (see (|25f) below). Then ljl2|) becomes 

oo a& 



r _ _ _ _ /*°° 

/ f(x)d,, = J2Y, (r). (18) 

•^ K " fe=0 ;=i ^° 



The main idea is simple and consists in replacing in formula Ijl8(l the non- 
negative univariate measures /ij! ; by their univariate GauB- Jacobi quadratures 5 



of order 2s — 1. Let ?/> 1 be the inverse map of rp and put ajfj — f^fe 8 ;) 
Then we obtain a pseudo-positive definite functional by setting 

oo gfc ,.00 

(/) : = E E / /w (0 '•"^S to ( 19 ) 

fc=0 i=l 
oo afc ^.oo 

= EE/ /m(v*)*-^Sw. 



fc=0 1 = 1 



As we have made it clear above, for / e C [x\,x%, we have (r) r - * = 

Pfc,z ( r2 ) 7 where pj.^ are polynomials. Hence, the functional T' s ) is well-defined 
on C [xi,X2, ■■■,Xd] , since then the series is finite. In fact, the most important 
thing is to find a condition on the measure fi which provides convergence of the 
series in (|19|1 for the class of continuous, polynomially bounded functions /. 

Since this is a very central result of our paper, let us give the main argument 
for proving the convergence of the series in (|19f) in the important case when 
all measures ; have their supports in the compact interval [0, R] . For the 
Laplace- Fourier coefficient, defined in JSJ, we have the simple estimate 

\fk,l {r)\<C max \f (x)\ for < r < R, 

\x\<R 

based on the Cauchy inequality and the orthonormality of {Yfc,z (&)} . Hence, 



h.i {r)r- k da\:Ur) 



<Cmax|/(x)| / r- k da { °](r) 

\x\<R ' 



and 



T« (/) < C max 1/ (*)l E E / r ~^S 

\x\<R * — ' ' — ' in 



(r). (20) 

k=0 1 = 1 

Now here is the crux of the whole matter: for the convergence in (|19fl it would 
suffice to prove the inequality 

poo />oo 

r- k d*l s ] (r) < / r- k dfi k j (r) . (21) 
Jo 

The famous Chebyshev extremal property 6 of the Gaufi-Jacobi quadrature pro- 
vides us with a proof of (|21l) . So we see that the convergence of the series in 

5 The upper index s will indicate the cardinality of the support. 

6 This has been proved by A. Markov EH1 ar >d T - Stieltjes, cf. Chapter 4] and I2H 

Chapter 3]. 
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I|19|l is a consequence of the summability condition (|llfl with N = 0. Further 
note that (|20l) shows that is a continuous functional: by the Riesz repre- 
sentation theorem we infer the existence of a signed measure with support 
in the closed ball B R := {x £ R n : \x\ < R} such that 

T^(f) = f f(x)da^(x) 

JB R 

for all continuous functions / : Bji — > C. Moreover, the component measures of 
the pseudo-positive measure are exactly the univariate measures The 
precise result is contained in Theorem [201 in Section In the case when not all 
measures fik,i have their supports in a compact interval [0, R] the argumentation 
has to be modified, and one needs (fill) for all N > 0. The details are provided 
in Section |21 and Section 0] 

The exactness of the Gaufi- Jacobi quadratures vf) for polynomials of degree 
< 2s — 1 implies that and /i coincide on the set of all polynomials P such 
that A 2s P — 0. This is due to the fact that in the Laplace-Fourier expansion 
Q the coefficients are given by fkj (r) = r k pk,i (r 2 ) where pt,i are polynomials 
of degree 2s — 1. For that reason we call the measure the polyharmonic 
Gaufi- Jacobi measure or the polyharmonic Gaufi- Jacobi cubature of order s. 

In the following we want to discuss the properties of the polyharmonic 
Gaufi-Jacobi cubature and it is natural to compare them with those of the 
univariate Gaufi-Jacobi quadrature. Among the various existing quadratures 
(e.g. Newton-Cotes quadratures), the Gaufi-Jacobi quadrature has the eminent 
property that the weights are positive. This in turn is the key to prove the 
convergence of the quadrature (see e.g. the discussion in ^3 p. 353] based on 
the theorems of Polya and Steklov). 

Property 1 (Stieltjes) For every continuous function f the Gaufi-Jacobi quadra- 
ture fdv^ converges to fdv, when s tends to infinity. 

A second important property of the Gaufi-Jacobi quadrature is the error 
estimate due to A. Markov (see ^3 p. 344]) 

Property 2 (Markov) Let v be a non-negative measure on [a,b] whose support 
has cardinality > s. Then for any 2s-times continuously differentiable function 
f : [a, 6] — > M there exists £ 6 (a, b) such that 

f f («) dv (t) - J" f (t) di/M (t) = J^yf (2s) (0 J \Qs (t)\ 2 (22) 

where Q s is the s-th orthogonal polynomial with respect to v, with leading coef- 
ficient 1. 

It is an amazing and non-trivial fact that properties ^) and [21 have analogs 
for the polyharmonic Gaufi-Jacobi cubature although the approximation mea- 
sures cr^ are in general signed measures. In Theorem 1221 we show that C s (/) 
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converges to / fd/x for every continuous function / : R" — * C. This property 
implies the numerical stability of our cubature formula. In Section [5] we prove 
an estimate for the difference 

M (/)-TW (/) 

for functions / G C 2s (Mr) by their derivatives in the ball Br based on Markov's 
error estimate. 

Let us outline the structure of the paper: In Section [21 we introduce the no- 
tions of pseudo-positive definite functional and pseudo-positive measure, and 
we prove basic results about them. In Section |3| we consider the determinacy 
question. In Section 0] the polyharmonic Gaufi-Jacobi cubature formula is pre- 
sented in detail. Section [5] is devoted to a multivariate generalization of the 
Markov's error estimate for the polyharmonic Gaufi-Jacobi cubature. 

Section|H|and[3are devoted to definition of polyharmonic Gaufi-Jacobi cuba- 
tures in other domains with symmetries as the annulus and the cylinder (periodic 
strip) . In Section H3 we construct a Gaufi-Jacobi cubature for pseudo-positive 
measures with support in a closed annulus A p .r which is exact on the space of 
all functions continuous on the closed annulus A Pt u and polyharmonic of order 
2s in the interior. While the case of the annulus is somewhat similar to that 
of the ball, we have to introduce a new notion of pseudo-positivity in the case 
of the cylinder (periodic strip) in Section [7| in order to obtain cubatures which 
preserve polyharmonic functions of order 2s. Moreover it is not possible to use 
in the proof the usual univariate Gaufi-Jacobi quadratures; instead we need the 
existence of quadratures of GauB-Jacobi-type for Chebyshev systems (Theorem 
I33[) . The analog to the crucial inequality l|21|) follows from the Markov-Krein 
theory of extremal problems for the moment problem for Chebyshev systems. 7 

In Section |H1 we give explicit examples illustrating our results and provide 
miscellaneous properties of pseudo-positive measures. In the last Section El 
we discuss shortly aspects of numerical implementation and some background 
information about the polyharmonicity concept. 

Finally, let us introduce some notations: the space of all continuous complex- 
valued functions on a topological space X is denoted by C (X) . By C c (X) we 
denote the set of all / € C (X) having compact support. Further C po i (R d ) is 
the space of all polynomially bounded, continuous functions, so for each / 6 
Cp i (R rf ) there exists N E N , such that \f (x)\ < C N {l + \x\) N for some 
constant Cn (depending on / ) for all x 6 M d . Further, we define an useful 
space of test functions 

N a k 

C x (R d ) := CM) F M (x):N€ N and e C [0, oo)}. (23) 

k=Q 1=1 

which can be rephrased as the set of all continuous functions with a finite 
Laplace-Fourier series. Moreover we set 

C c x (E rf ) := C x (R d ) n C c (R d ) . (24) 

7 We use the name "Markov-Krein theory" following 24' Chapter 3], while in 1321 this is 
called " Chebyshev-Markov problem". 
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We need some terminology from measure theory: a signed measure on R d is 
a set function on the Borel er-algebra on R d which takes real values and is 
er-additive. For the standard terminology, as Radon measure, Borel er-algebra, 
etc., we refer to |Sj. By the Jordan decomposition p. 125], a signed measure 
[i is the difference of two non-negative finite measures, say [i = /i + — (JT with the 
property that there exist a Borel set A such that /i + (A) — and (R™ \ A) = 
0. The variation of (i is defined as := fi + + fi~ . The signed measure [i is 
called moment measure if all polynomials are integrable with respect to /i + and 
which is equivalent to integrability with respect to the total variation. The 
support of a non-negative measure /i on R rf is defined as the complement of the 
largest open set U such that /i (U) = 0. In particular, the support of the zero 
measure is the empty set. The support of a signed measure a is defined as the 
support of the total variation |<r| = ct + + cr_ (see p. 226]). Recall that in 
general, the supports of a+ and er_ are not disjoint (cf. exercise 2 in ^] p. 
231]). For a surjective measurable mapping ip : X — > Y and a measure v on X 
the image measure v v on Y is defined by 

^(B):=y{ V - 1 B) (25) 

for all Borel subsets B of Y. The equality f x g (<p (x)) dv (x) = J Y 9 (y) dv 1 * (y) 
holds for all integrable functions g. We use the notation 0Jd-i for the surface 
area of the unit sphere, so 

Wd-i ■= f IdO. (26) 

By B R wc denote the closed ball {i £ 1" : |x| < R} . For < p < B < oo wc 
define the closed annulus by 

A p ,r := {x € R d : p < \x\ < R} . (27) 

For a subset B of R" the interior is denoted by B° . Further, we denote the 
closed and the open interval respectively by [a, b] = {x € R : a < x < b} and 
(a,l)) = {iel:s<i<i)}. 

2 The moment problem for pseudo-positive def- 
inite sequences 

Let T : C [x\, X2, — ► C be a linear functional. For any solid harmonic 
polynomial Ykj (x) we define the component functional T^j by 

T k ,i (p) := T c (p(\x\ 2 )Y Kl (x)) for every p £ C [a*] . (28) 

Let us give the precise definition of pseudo-positive definiteness, which we al- 
ready mentioned in the introduction: 
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Definition 3 A sequence c — {c a } a ^d , or the associated functional T Cl is 
pseudo-positive definite if for every k £ No and I = l,...,<ifc the sequences 
{ c j' fe '^}ieN an d { c ^+i}j'eNo defined in are positive definite. Clearly this is 
the same to say that T ki i (p* (t)p(t)) > and T k ,i (t ■ p* (t)p(t)) > for every 
P (t)eC[ Xl }. 

First we recall the following result which may be found e.g. in 6 or 44 . 

Proposition 4 The Laplace- Fourier coefficient fk,i of a polynomial f given by 
(5)) is of the form fk,i if) = r k pk.i (r 2 ) where pk,i is a univariate polynomial. 
Hence, the Laplace- Fourier series ^ is equal to 

/(^^ Py (|x| 2 )r w (x). (29) 

k=0 1=1 

Equality (|29|l is a reformulation of the Gaufi decomposition of a polynomial 
which we have provided in (J3J). 

The next two Propositions characterize pseudo-positive definite sequences: 

Proposition 5 Let c = {c a } a6N d be a pseudo-positive definite sequence and 
T c its associated functional. Then for each k £ No, I = l,...,dk, there exist 
non-negative measures o~k,i with support in [0, oo) such that 

dcg/ a fc „oo 

r c (/)=EE| fkAr)r- k da ktl (r) (30) 
fc=0 i=i •'° 

holds for all f 6 C [x\, X2, x^] where /jy (r), k 6 No, I = 1, ...,dfc, are the 
Laplace- Fourier coefficients of f. 

Proof. By the definition of pseudo-positive defmiteness, T k ,i (p* (t)p(t)) > 
and Tk t i(t ■ p* (t)p(t)) > for each univariate polynomial p (t) where the 
component functional T k ,i is defined in (|28|l . By the solution of the Stieltjes 
moment problem there exists a non-negative measure /j, k j with support in [0, oo) 
representing the functional Tfc,;, i.e. satisfying 

/■oo 

T k ,l(p)= / p(t)dfi k ,l{t) for every p e C [t] . (31) 
Jo 

Let now <p : [0, oo) — > [0, oo) be defined by 99 (t) = \ft. Then we put cr k j '■= (J%, 
where ^ ; is the image measure defined in l(2"B|) . We obtain 

h(t)dn k>l (t) = / h(r 2 )d^ kl {r). (32) 
Jo 

Now use (|29|) . the linearity of T and the definition of T k> i in i|28[l . and the 
equations (|31|1 and l)32|l to obtain 

deg/ a fe dcg/ a k .qq 

T c (/)=££ r <M (?*.o - E E / (^ 2 ) w • 

fc=0 1=1 fc=0 {=1 ^° 
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Since p k j {r 2 ) = r~ h f kt i (r) the claim (I3(J[1 follows from the last equation, which 
ends the proof. ■ 

The next result shows that the converse of Proposition's also true; not less 
important, it is a natural way of defining pseudo-positive definite sequences. 

Proposition 6 Let o~ kj i, k £ No , I = 1 , . . . , , be non-negative moment measures 
with support in [0, oo) . Then the functional T : C [xx,X2, Xd] — > C defined by 

dcg/ a h roo 

fk,l (r) r- k da ktl (33) 



T (f) -EE/ 

k=o i=i Jo 



is pseudo-positive definite, where fk,i (r), k 6 No, I = 1, ■■■,a k , are the Laplace- 
Fourier coefficients of f. 

Proof. Let us compute T k j (p) where p is a univariate polynomial: by defi- 
nition, Tk.i (p) — T (p(\x\ 2 )Yk.i (x)\. The Laplace-Fourier series of the function 
x I— > |as| ■'pdx] )Y k j (x) is equal to r 2j p (r 2 ) r^Ife,; (0), hence 

T k>l (t?p{t)) =T(\x\ 23 p{\x\ 2 )Y Kl {x)) = I™ rip{r 2 )da k j 



for every natural number j. Taking j • = and j = 1 one concludes that 
T k ,i (p* (t) p (t)) > and T k j (tp* (t) p (t)) > for all univariate polynomials 
p, hence T is pseudo-positive definite. ■ 

The pseudo-positive definiteness is defined for afunctional on C [xi, X2, Xd] ■ 
Now we introduce the concept of pseudo-positivity of a measure: 

Definition 7 ^4 signed measure [i on R" is called pseudo-positive if 

h(\x\)Y k<l (x)dfi(x)>0 (34) 



holds for every non-negative continuous function h : [0, oo) — * [0, oo) with com- 
pact support. 

At first we need some basic properties of pseudo-positive measures. 

Proposition 8 Let p, be a pseudo-positive moment measure on R d . Then there 
exist unique moment measures fj, k j defined on [0, oo) , which we call component 
measures, such that 



h(t)da ktl (t)= / h(\x\)Y K i(x)du (35) 
holds for all h E C po i [0, oo). Further for each f e C p x , (R d ) 

/ /(x)d/i = VY) / f k ,i(r)r- k du kjl . 



fc=0 z=i 
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Proof. By definition of pseudo-positivity, M k ,i (h) := / ffid h (\x\) Y k j (x) d\i 
defines a positive functional on C c ([0,oo)). By the Riesz representation the- 
orem there exists a unique non-negative measure fi k _i such that M/.J (h) — 
J °° h (t) dfik,i for all h £ C c ([0, oo)) . We want to show that l|3l)|) holds for all 
h £ Cp i [0, oo). For this, let ur : [0, oo) — > [0, 1] be a cut-off function, so ur is 
continuous and decreasing such that 

u R (r) = 1 for all < r < R and u R (r) = for all r > R+ 1. (36) 

Let /i £ Cp j [0, oo) . Then it^/i G C c ([0, oo)) and 

Ur (*) h (t) dnk,i = / UK(|a:|)M|a:|)y*,j(aOdft. ( 37 ) 

Note that (t) /i (t)\ < \ur+i (t) h (t)\ for all t £ [0, oo) . Hence by the mono- 
tone convergence theorem 



\h(t)\dii k ,i = Jim / |u R (i)/i(i)|c^ fe ,i. (38) 
On the other hand, it is obvious that 



/ 



u R {\x\)\h(\x\)\Y k j (x)d f i 



< / \h{\x\)Y k}l (x)\d\fi\. (39) 



The last expression is finite since /i is a moment measure. From (JSSj, (|37H 
applied to \h\ and (|3^|l it follows that \h\ is integrable for /j, k j. Using Lebesgue's 
convergence theorem for [i and (|37|l it is easy to that (|35|l holds. For the last 
statement note that each / £ C* ol (R d ) has a finite Laplace-Fourier series, and 
it is easy to see that the Laplace-Fourier coefficients f k ,i are in C po i [0, oo), see 
|g2J below. ■ 

The next theorem is the main result of this section and it provides a simple 
sufficient condition for the pseudo-positive definite functional on C [xi, xi, Xd] 
defined in (|33|l to possess a pseudo-positive representing measure. Let us note 
that not every pseudo-positive definite functional has a pseudo-positive repre- 
senting measure, see Section [3] for an example. 

Theorem 9 Let a kl i, k £ Nq, I = 1, ■ ■■,a k , be non-negative measures with sup- 
port in [0, oo) such that for any N £ Nq 

Cn:=Y.Y. rN r- k d<J kA < oc . (40) 
fc=0 1 = 1 ^° 

Then for the functional T : C \x\,x<i, ■■■,Xd] —* C defined by there exists a 
pseudo-positive, signed moment measure a such that 

T(f) = J fda for all f £ C[x 1 ,x 2 , ...,x d ] ■ 
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Remark 10 1. If the measures a k ,i have supports in the compact interval [p,R] 
for all k £ No, I = 1, <Xfc, then the measure a in Theorem\Qhas support in the 
annulus {x £ M. d : p < \x\ < i?} . 

2. In the case of R < oo , it obviously suffices to assume that Co < oo 
instead of Cm < oo for all N £ No . 

3. The proof of Theorem^ shows that o~k,l * s equal to the component measure 
induced by a with respect to the solid harmonic Y k) i (x) . 

Proof. 1. We show at first that T can be extended to a linear functional T 
defined on C po i (R d ) by the formula 

OO Qfc «oO 

T (/) : =EE/ hAr)r~ k da k , (41) 
k=a i=i 

for / £ Cp i (R d ) , where f k ,i (r) are the Laplace-Fourier coefficients of /. In- 
deed, since / £ C po i (R d ) is of polynomial growth there exists C > and TV £ N 
such that |/ (a:) | < (7(1 + \x\ N ). If follows from © that 



\f k ,l(r)\<C(l + r N )^M-jJ \Y ktl {e)\ 2 de = C{l + r N )^U, (42) 

where we used the Cauchy-Schwarz inequality and the fact that Y k ,i is orthonor- 
mal. Hence, 

n oo 

\f k ,l {r)\r- k da kI < JUZHC \ (l + r N ) r~ k da k ,i. 



By assumption l|4U|) the latter integral exists, so f k ,i (r) r k is integrable with 
respect to a k j. By summing over all k, I we obtain by (|4U|I that 



EE 

fc=0 1=1 



f k ,i (r) r fc dcr M 



o 



< 00, 



which implies the convergence of the series in <|41ll . It follows that T is well- 
defined. 

2. Let To be the restriction of the functional T to the space C c (R d ) . We will 
show that To is continuous. Let / £ C c (R d ) and suppose that / has support in 
the annulus {x £ R d : p < \x\ < i?} (for the case p = this is a ball). Then by 
a similar technique as above 

\fk,l (r)\ < V^d-i max \f(x)\. 

p<\x\<R 

Using (|41|l one arrives at 

oo a k .R 

\T (f)\< max j/(i)|v^EE / r ~ k do~ ki i ■ (43) 
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3. First consider the case that all measures a^j have supports in the interval 
[p,R] with R < oo (cf. Remark HU|l . Then and the Riesz representation 
theorem for compact spaces yield a representing measure with support in the 
annulus {a; £ R d : p < \x\ < R} . The pseudo-positivity of p will be proved in 
item 5.) below. 

In the general case, we apply the Riesz representation theorem given in |SJ 
p. 41, Theorem 2.5]: there exists a unique signed measure a such that 

To (<?) = / Qda for all 9 £ C c (R d ) . 

JR d 

4. Next we will show that the polynomials are integrable with respect to 
the variation of the representation measure a. Let a = a+ — o- be the Jordan 
decomposition of u. Following the techniques of Theorem 2.4 and Theorem 2.5 
in p. 42], we have the equality 

/ g (x) da + = sup {T (h) : h G C c (R d ) with < h < g} (44) 

which holds for any non-negative function g E C c . Let ur be the cut-off 
function defined in (|36|) . We want to estimate J Rd g (x) da+ for the function g := 

\x\ N u R {\x\ 2 ). In view of g2J|, let h e C c (R d ) with < h (x) < \x\ N u R (\x\ 2 ) 
for all x £ M d . Then for the Laplace-Fourier coefficient hk,i of h we have the 
estimate 

\h k ,i(r)\<J [ \h{r6)\ 2 d9 x f \Y k<l (6)\ 2 dO < r N u R (r 2 ) > /E£T. 
According to lUTl) 

n (h) < \t (h)\ < v^^Y, Yl / rN ^ kd<j k,i =■■ d n . 

k=0 1=1 •'° 

From fLl|) it follows that J Rd \x\ N u R (\x\ 2 )da+ < Dn for all R > (note that 
Dn does not depend on R ). By the monotone convergence theorem (note that 
ur (x) < Uft+i (x) for all x £ M. d ) we obtain 

/ 1x1^^0"+= lim / \x\ N U]i(\x\ 2 )da+ < Dn- 

Similarly one shows that J Rd \x\ N da^ < 00 by considering the functional S = 
—Tq. It follows that all polynomials are integrable with respect to <r+ and <r_. 
Using similar arguments it is not difficult to see that for all g £ C x (R d ) n 
Cp i (R d ) 

f g(x)da = T(g). (45) 
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5. It remains to prove that a is pseudo-positive as given by Definition 
Let h £ C c ([0, oo)) be a non-negative function. The Laplace-Fourier coefficients 
fk',l> of / (x) ■= h (\x\) Y kjl (x) are given by } k - v (r) = 8 kk ,S w h{r) r k and by gS} 
it follows that 



h(\x\)Y kJ (x)da = T(f) = fk,i(r)r-"dtT k ,i= h(r)da kil . 

Jo Jo 

Since a k ,i are non-negative measures, the last term is non-negative. According 
to definition J3IJ, a is pseudo-positive. The proof is complete. ■ 

The following is a solution to the modified moment problem as explained in 
the introduction. It is an immediate consequence of Theorem [5] 

Corollary 11 Let T : C [x\, X2, x<j] — ► C be a pseudo-positive definite func- 
tional. Let a kt i, k € No, I — 1, <2fc, 6e i/ie non-negative measures with supports 
in [0, oo) representing the functional T as obtained in Proposition^ If for any 
NeN 



oc 



r N r- k da ktl < oo, (46) 

fe=0 1=1 ^° 

i/ien i/iere exists a pseudo-positive, signed moment measure a such that 
T(f) = J fda for all f e C[xi,X2, ■■■,Xd] ■ 

It would be interesting to see whether the summability condition (|46J) may 
be weakened, cf. also the discussion at the end of Section [H] 

By the uniqueness of the representing measure in the Riesz representation 
theorem for compact spaces we conclude from Theorem 

Corollary 12 Let fi be a signed measure with compact support. Then fi is 
pseudo-positive if and only if fj, is pseudo-positive definite as a functional on 

C[x 1 ,X2,...,Xd\. 

Let us remark that Corollary ^| does not hold without the compactness 
assumption which follows from well known arguments in the univariate case: 
Indeed, let v\ be a non-negative moment measure on [0,oo) which is not de- 
termined in the sense of Stieltjes; hence there exists a non-negative moment 
measure vi on [0, oo) such that v\ (p) = Vi (p) for all univariate polynomials. 
Since v\ ^ V2 there exists a continuous function h : [0, oo) — > [0, oo) with com- 
pact support that v\ (h) ^ v-i (h) . Without loss of generality assume that 

/•OO poo 

/ h (r) dv t - / h (r) dv 2 < 0. (47) 
Jo Jo 

For i = 1,2 define fXi — dOdvi, so for any f £ C (R d ) of polynomial growth 



//•OC P 
fdm= / / f{r6)d6d Vi 
JO J§d-1 
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For a polynomial / let fo be the first Laplace-Fourier coefficient. Then J fdfii — 
Jo fo ( r ) dvi for i = 1,2. Since vi (p) — v% (p) for all univariate polynomials it 
follows that J fdfii — J fd{i 2 for all polynomials. Then fi := /ii — /i 2 is a signed 
measure which is pseudo-positive definite since [i (P) = for all polynomials P. 
It is not pseudo-positive since /j,q (h) = J h (\x\) d/i < by 1|47|> . 



3 Determinacy for pseudo-positive definite func- 
tionals 

Let M* (R d ) be the set of all signed moment measures, and M£ (R d ) be the set 
of non-negative moment measures on R d . On M* (R d ) we define an equivalence 
relation: we say that a ~ /i for two elements a,[i E M* (R rf ) if and only if 
J Rd /da = J Rd fd/J, for all / € C [x\, x 2 , x d ] . 

Definition 13 Let fi E M* (R d ) be a pseudo-positive measure. We define 

= {cr E M* (M. ) : a is pseudo-positive and a ~ //} . 

VFe say i/iai i/ie measure [i E M* (R d ) is determined in the class of pseudo- 
positive measures if V M /las on/y one element, i.e. is equal to {/i} . 

Recall that a positive definite functional tf> : V\ — * R is determined in the 
sense of Stieltjes if the set 

J^f 1 := |r E M* ([0, oo)) : jf r m dr = (r m ) for all m G N | (48) 

has exactly one element, cf. 9, p. 210]. 

According to Proposition^ we can associate to a pseudo-positive measure /z 
the sequence of non-negative component measures fik,i, k G No, / = 1, •-, a k with 
support in [0, oo) . The measures fj, k ,i contain all information about fx. Indeed, 
we prove 

Proposition 14 Let /i and a be pseudo-positive measures and let ^ ; and o~k,i 
be as in Proposition^ If fJ-k,i = &k,i for all k E No, I = 1, .., a k then ji = a. 

Proof. Let h E C c [0, oo) . Then, using the assumption fi kt i = Uk,i, we obtain 
h(\x\)Y kt i(x)dfj,= / h(t)dfjtk,i= / h(\x\)Y kil (x) da. 



Since each / G C c x (R d ) is a finite linear combination of functions of the type 
h (\x\) Yk,i (x), we obtain that / Rd fd\i = J Rd fda for all / G C c x (R d ) . We apply 
Proposition El to see that /i is equal to a. ■ 

The following result is proved in [§1 Proposition 3.1]: 
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Proposition 15 Let \l and a be signed measures on M. d . If J Rd fd[i — JL d fda 
for all f € C c x (R rf ) , then fi is equal to a. 

We can characterize in the case that only finitely many are nonzero. 

Theorem 16 Let fi be a pseudo-positive measure on 1™ such that fik.i — for 
all k > k 0l I = 1, dfc. Then Vu is affinely isomorphic to the set 



poo 
JO 



c c/p M < oo} (49) 



where the isomorphism is given by a i — ► (^fci) ^ e m «P 

^ : [0, oo) — > [0, oo) is defined by ip (t) = t 2 , cf. |SSp . 

Proof. Let cr be in V^. Let cr^; and jikj be the unique moment measures 
obtained in Proposition [S] Then 

h (t) daft = [ h (t 2 ) da k ,i = f h{\x\ 2 )Y kA (x) da (x) 

JO JR n 

for all h £ Cp i [0, oo) , and an analog equation is valid for fjL/,,1 and fx. Taking 
polynomials h (t) we see that ak,l & W s * 1 using the assumption that // ~ cr. 

Using a simple approximation argument it is easy to see from (|35|) that 



Since x i — > Yfe } / ( t^t ) is bounded on E n , say by M, we obtain the estimate 



Jo 



< M / Idler! < oo. 



It follows that (erf A is contained in the set on the right hand 

side in (gSJ. 

Let now € W% % be given such that L t~^ k dpk.i < oo for k = 1, .., ko, I = 

1, ...,afc. Define cr^ j = and cr^; = for fc > &o. Then by Theorem [5] there 
exists a measure reVj, such that Tfe^ = cr^;. This shows the surjectivity of the 
map. Let now a and r are in with erf, = rt, for k = 1, ..,/co,/ = 1, ...,afc. 
The property cr € implies that cr^ ; G for all fc 6 No,/ = 1, ...,cifc, 

hence cr^ , = for k > ko, and similarly rt, = 0. Hence (Jk,l — T k.i for ah 
€ No, I = 1, dfc, and this implies that a — t by Proposition ■ 
The following is a sufficient condition for a functional T to be determined in 
the class of pseudo-positive measures. 



18 



Theorem 17 Let T : C [x%, xi, ■ Xd] — > R be a pseudo-positive definite func- 
tional. If the functionals Tkj : C [xi] — > C are determined in the sense of 
Stieltjes then there exists at most one pseudo-positive, signed moment measure 
yU on R d with 

T(f)= [ fdfi for all f £C[x u x 2 ,...,x d ]. (50) 

Proof. Let us suppose that fi and a are pseudo-positive, signed moment 
•as 
that 



measures on R d representing T. Taking / = |x| 2 Yjy (x) we obtain from |5 



\x\ 2N Y Ki (x) dp, = T Kl (t N ) = / \x\ 2N Y Kl (x) da 



for all N £ No- Let fik,i and &k,l as in Proposition 03 and consider if> : [0, oo) — > 
[0,oo) defined by ip (t) = t 2 . Then the image measures /ij^ and are non- 
negative measures with supports on [0, oo) such that J °° t N dpf l — Tk,i (t N ) = 

fo° t N dcrf j. Our assumption implies that = af [, so = &k,l- Proposition 
E| implies that /i is equal to a. ■ 

In the following we want to prove the converse of the last theorem, which 
is more subtle. We need now some special results about Nevanlinna extremal 
measures. Let us introduce the following notation: for a non-negative measure 
<P £ Ml (R) we put 8 

[4>] := {a £ M* (R) : a ~ <f>} . 

Proposition 18 Let v be a non-negative moment measure on R with support in 
[0, oo) which is not determined in the sense of Stieltjes, or applying the notation 
W„ tl ^ {y} . Then there exist uncountably many a £ W„ t% such that 
r k da < oo for all k £ Nq. 

Proof. In the proof we will borrow some arguments about the Stieltjes 
problem as given in or ES ■ As in the proof of Proposition 4.1 in |59"| 
let </? : (—00,00) — > [0,oo) be defined by y (x) = x 2 . If A is a measure on R 
define a measure A~ by A~ (A) := A (—A) for each Borel set A where —A := 
{—x : x £ A} . The measure is symmetric if A - = A. For each r £ W^ u define 

a measure f := | [r v + (t v ) ^ which is clearly symmetric, in particular is 

symmetric. As pointed out in the map r : W 7 ^* 1 — > \v\ is injective and the 
image is exactly the set of all symmetric measures in the set \y] . The inverse 
map of ~ defined on the image space is just the map a — > a v . 

It follows that v is not determined, so we can make use of the Nevanlinna 
theory for the indeterminate measure v, see p. 54 in pQ. We know by for- 
mula II.4.2 (9) and II. 4. 2 (10) in £Q that for every f e 1 there exists a unique 
Nevanlinna-extremal measure at such that 

da t (w) _ _A(z)t-C (z) 
, u- z ~ ~ B(z)t- D(z)' 

8 Here in order to avoid mixing of the notations, we retain the notation [<f>] from the one- 
dimensional case in UJ. 
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where A (z) , B (z) , C (z) ,D (z) are entire functions. Since the support of a t is 
the zero-set of the entire function B (z)t — D (z) it follows that the measure at 
has no mass in for t 7^ 0, and now it is clear that £7t([— S, 5]) = for t 7^ and 
suitable 8 > (this fact is pointed out at least in the reference 9, p. 210]). It 
follows that 

/>oo 

|u| da t < 00 (51) 

— 00 

since the function u 1 — ► |u| k is bounded on R \ [— 5, S] for each 6 > 0. Using 
the fact that the functions A (z) and B (z) of the Nevanlinna matrix are odd, 
while the functions B (z) and C (z) are even, one derives that the measure pt := 
^at + \a~ t is symmetric. Further from the equation A (z) D{z) — B (z) C (z) = 1 
it follows that pt =/= p s for positive numbers t 5^ s. By the above we know that 
Pt Pt- This finishes the proof. ■ 

Theorem 19 Let p be a pseudo-positive signed measure on M. d such that the 
summability assumption holds. Then contains exactly one element if 
and only if each pt , is determined in the sense of Stieltjes. 

Proof. Let pk,i be the component measures as defined in Proposition |SJ 
Assume that — {p} but that some r := pf g lg is not determined in the sense 
of Stieltjes where ip (t) — i 2 for t G [0, 00) . By Proposition ^] there exists a 
measure a £ W^ tl such a 7^ r and J °° r~ k da < 00. By Theorem [5] there exists 
a pseudo-positive moment measure p representing the functional 

00 a k />oo y* 00 

T (/) : = E E / hAr)r- k dp k ,+ f ko , lo {r) r~ k da*- 1 . 

k=a,k^k l=l,tyl 

Then p is different from p since a^ 7^ Pk ,i an d p € since a G 

This contradiction shows that pf z is determined in the sense of Stieltjes. The 

sufficiency follows from Theorenili^l The proof is complete. ■ 



4 Polyharmonic Gaufi- Jacobi cubatures 

In this section we will prove the main result of the paper, the existence of the 
polyharmonic Gaufi- Jacobi cubature of order s. The proof is based on applica- 
tion of the famous Chebyshev extremal property of the Gaufi- Jacobi measure. 

Theorem 20 Let < p < R < 00. Let p be a pseudo-positive signed measure 
with support in the closed annulus A p ^ such that 

SSImh)^ 00 ' (52) 

Then for each natural number s there exists a unique pseudo-positive, signed 
measure a^ with support in A Pt u such that 
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(i) The support of each component measure a k ] of (defined by \35\) ) 
has cardinality < s. 

(ii) f Pdfi = J Pda^ for all polynomials P with A 2s P = 0. 

Proof. By Proposition [S] the following identity holds 

/. OO dk />OC 

/ f(x)d^(x)=J2J2 h,l{r)r- k d m {r) (53) 
Jm" k=0 l=x Jo 

for any / G C* ol (K d ) where fik,t (h) = J h (|x|) Ykj (x) dfi (x) . It is clear that 
fik,l has support in the interval [p, R) . If the cardinality of the support of /j,k,i 
is < s we define a k [ := fik,l- If the cardinality is strictly larger than s we define 

(s) 

a k [ as the non-negative measure such that 

R rR 

r 23 da[ s j (r) - / r 2 ^ dfi kJ (r) (54) 

for all j = 0, 2s— 1. The existence of a k s j is proved as follows: Let ip : [p, R] — > 
[p 2 , R 2 ~\ be the map ip (t) — t 2 . Then the image measure /if ; is a measure on 
[p 2 , R 2 ] and its support has clearly cardinality > s. Let fffj be the Gaufi-Jacobi 
quadrature of [i k { . From the Gaufi-Jacobi quadrature formula (|To)l (see also the 
footnote after it) follows that v k s J has support in the open interval (p 2 , R 2 ) and 

r R 2 



= \ tH^ 

J p 2 

for j — 0, 2s — 1. Now it is easily seen that a k s j := (^ii*) satisfies (|54H . 
We now define a functional T : C [xi, X2, Xd] —> C by putting 



deg/ a fc „jj 

T (s Hf) = EE/ /mW'-HI- ( 55 ) 

fe=0 2 = 1 "'P 



Let us show that 

P(x)dfi(x)=T^(P) (56) 



for all polynomials P with A 2s P (x) = 0. Indeed, according to l|29|l the Laplace- 
Fourier series of a polyharmonic polynomial of order 2s can be written as 

deg / afe 
fc=0 2 = 1 

and the univariate polynomials p k ,i (t) have degree < 2s — 1 (see e.g. |3 II 
Theorem 10.42, Remark 10.43]). Combining (JSTJ) with J55J), {32) and {53J gives 
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Now we want to prove that can be represented by a signed measure. 
We claim that 

R pR 

r- k da$ (r) < / r- k d m (r) < oo (58) 
p Jp 

for all k £ No, I = 1, ...,cifc. If cr^] = ^ fcj ; there is nothing to prove, so we can 
assume that the cardinality of the support of /ik,i is bigger than s. After taking 
the image measures under the map ijj we see that we have to prove 

*-**<&$(*)< f H*d/*J,(t). (59) 

For p > this follows from the Chebyshev extremal property of the GauB- 
Jacobi measure (see e.g. [33 Chapter 4, Theorem 1.1]) applied to the function 
/(*) := t-2 k . The same result works in the case p = but due to the singularity 

( s) 

of / we have to use essentially the fact that all points of v k [ are in the open 
interval (0,i? 2 ) and to apply Remark 1.2 in Chapter 4 of [22]. 9 We have only 
to check that the assumptions on / are satisfied: / is non-negative, 2s times 
differentiable on the open interval (p 2 ,R 2 ), and for t £ (p 2 ,i? 2 ) 

^ = '- 1 ) 2 ' ¥ + • + 2 » - > »■ <»> 

By our assumption (|52|) and by (|58() we can apply Theorem [!|] and (i) is 
proved. Property (ii) follows from (|56[) which we have proved above. 

Let us prove the uniqueness of a^ s \ Assume that r is a signed pseudo- 
positive measure with compact support, and with properties (i) and (ii). Since 
t is pseudo-positive there exists by Proposition |S| univariate measures Tk t i such 
that 

r h (t) dr k j = [ h(\x\)Y kl i (x) dr (x) (61) 

for any polynomially bounded continuous function h. Since A 2s (^\x\ 2 ^ Y k ,i (x) 
for j = 0, 2s — 2, we infer that 

\x\ 2 ' Y kA (x) dr (x) = J \x\ 2 ' Y kA (x) dfi (x) . 

Hence / t 2j dTk,i = J t 2j dpk,i for j — 0, 2s — 2, so 

t'drt, = ( Vd4,. (62) 



9 It is curious that Stioltjes proved that the Gaufi— Jacobi quadrature measure solves a three- 
dimensional spherically symmetric extremal problem with a singular function / (t) = see 

v t° 

the complete description in 1321 Chapter 4.2, formula (2.6) ]. 
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By property (i) the support r^; has cardinality < s, hence r kl has cardinality 
< s. The uniqueness of the Gaufi-Jacobi quadrature shows that r kl is equal to 

(s) 

v k [ which means that Tk,i = &k,l- If the support of [ik,i has less than s points 
then crjj. s ] in our construction is defined to be fik,l- From l|62[l one can derive 



that r kl has the same support set as \x kl and finally that Tk lt 
Proposition E| yields r = ij( s \ ■ 



A*fc,z 



'k,i ■ 



Definition 21 The measure constructed in Theorem \2(A will be called the 
polyharmonic Gaufi-Jacobi measure of order s for the measure /x. 

The following is an analog to the theorem of Stieltjes about the convergence 
of the univariate Gaufi-Jacobi quadrature formulas. 

Theorem 22 Let < R < oo and let be the polyharmonic Gaufi-Jacobi 
measure of order s for the measure /x, obtained in Theorem \2(A Then 



f (x) da^ — ► f (x) da for s — > oo 



holds for every function f € C (Bji). 

Proof. Item (ii) of Theorem [3U] implies that for any polynomial P the 
convergence T^ (P) -> P holds for s — > oo. Theorem 14.4.4 in [E] shows that 
the convergence T^ s > (/) — * / carries over to all continuous functions / : Br — > C 
provided there exists a constant C > such that for all natural numbers s and 

& llfeC(B R ) 

T«(/) <C max 1/(201. 

But that is just estimate (|43|l in the proof of Theorem[5] The proof is complete. 
■ 

Using the same techniques as in Theorem 1201 we may prove a generalization 
of the Chebyshev extremal property of the Gaufi-Jacobi quadrature: 

Theorem 23 Let < p < R < oo and let u be a pseudo-positive signed measure 
with support in A p r satisfying the summability condition $5fy) and let <j( s J be 
the polyharmonic Gaufi-Jacobi measure of order s. Let f £ C 2s (R d ) be such 
that for all t € (p 2 , R 2 ) holds 



d 2s 
dt 2 ^ 



fk.i (Vtj t- 1 ^ 



> 0, 



for all k £ No, I = 1, 2, a^. Then the following inequality 

J f(x)dtrW < J f(x)du 

holds. 
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Although the measures a k s [ are based on point evaluations, it is clear that 
our approximation measures <jW are not point evaluations. The following gives 
a description of the support of the polyharmonic Gaufi-Jacobi measure <jW 
when only finitely many measures fik,i are non-zero: 

Proposition 24 Suppose that o~k,h k = 1, ...,fco, I = 1, ■■■,ak, are non-negative 
measures with finite support in the open interval (0, oo) and suppose that o~kj — 
for all k > ko ■ Then the support of the representing measure of the functional 
T : C [xi, X2, ■ Xd] —> C defined by is contained in the union of finitely 

many spheres with positive radius and with center 0. 

Proof. Let Sk,i be the finite support of cr^l an d let S be the union of all 
sets S k ,i with k = 1, ko, I = 1, a k . Let S — {x e R d : |x| S 5} . We show 
that the support of a is contained in S. Indeed, one can estimate as in 

fcn ak 

\T{P)\ <max|F(x)|V^T^X] / r' k da ktl (t) . 

By the Riesz representation theorem, T can be represented by a signed measure 
a which has support in the compact set 5*. ■ 

Remark 25 An interesting characteristic feature of the classical Gaufi-Jacobi 
quadrature measure is the minimality of its support among all non-negative 
measures which are exact of the same degree 2s — 1. One might see above some 
analogy with this phenomenon if one considers fi k ,i (i") and <ri j (r) as measures 
defined on the space K = {(k, I, r) : k 6 No, I = 1, 2, cjfc, r £ [0, oo)} . 



5 Markov type error estimates 



In this section we want to give an error estimate for our cubature formula. The 
proof is based on the Markov estimate for the Gaufi-Jacobi measure provided 
in 1(52)1, 

Let s £ NU{oo} . For an open subset U of M. d we denote by C s (U) the space 
of all / G G (U) which are continuously differentiable in U up to the order s. 

Theorem 26 Let < p < R < oo and let ip : [0, oo) — > [0, oo) be defined by 
ip (t) = t 2 . Let be a pseudo-positive signed measure with support in A p> r sat- 
isfying the summability condition and let crW be the polyharmonic Gaufi- 
Jacobi measure of order s. Define for every f £ C (A p j{) the error functional 

E s (f) := / /(x)d/x(x)- / /(x)dcrW(x). 



If f E C 2s (A° R ) n C {A P . R ) then E s (f) is lower equal to 



p 

oo a k 



(2s) 



EE, 



sup 



;=0 1=1 p 2 <«<« 2 



dt 2s 



\QL (t)\ d4 
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Here Q s k , (t) is the orthogonal polynomial of degree s with respect to the measure 

fj,t,, normalized so that the leading coefficient is equal to 1; if the support of p k ,i 
has less than s points, Qf, , is defined to be 0. 

Proof. Since / £ C 2s (A° R ) n C (A Pt ft) it is easy to see that the Laplace- 
Fourier coefficients f k ,i € C 2s (p,R) H C[p,R]. Let pj k j ana - °~k,U k € No, I = 
1, afe, and er( s ) be as in Theorem 1201 From the definitions it follows 



oo a k . 



/k,j (r) r- fc ^ fc>i - / /*,, (r) r" 



fc,i ■ 



Further /jy (r) r fc is integrable with respect to fif.,1 since ff.j is continuous 
on [p, R] and condition (|52|) holds. Let us fix the pair of indices (k, I) . If the 
support of fik,i has less than s points we know that pL k ,i = o~ki- So assume 
that the support of (if.,1 has at least s points. Then the support of fij^ l has at 

least s points and in our construction vjf] is the Gaufi-Jacobi measure of p,f t . 
Consequently 

e (fk,l) ■■= f fk,l (r) r~ k dp Kl (r) - f f k<l (r) r^da^ (r) 
<-R 2 , , rR 2 



By the proof of Markov's error estimate (1221) given in (JHl) one easily obtains 
with g ky i (t) := f k ,i (Vt) t~? k the inequality 

e(f k j)<-^ sup 9^(0 f R \Qh(t)\ 2 d^ k .(t). 

(2S)\ p 2 <i<R 2 J p 2 

The proof is complete. ■ 

In the following we want to give a Markov type error estimates for holomor- 
phic functions /. We will need the following property which was observed in 

u 

Lemma 27 Let f € C°° (B° R ) . Then f kJ E C°° [0, R) and £^f kJ (0) = for 
m = 0, k — 1. 

Lemma 28 Let f be a holomorphic function on the open ball :— {w € C d : 
j \u>j\ 2 < t 2 } for t > 0. Let f k> i be the Laplace- Fourier coefficient of f 
given by J3|) and let p k .i (t) be defined by the equation f kt i (r) — p k ,i (r 2 ) r k for 
< r < t. Let p and t satisfy < t 2 < p < r. Then 



d s 

(*) 



hold for all s = 0, 1, 2, . 



2 — k | 

<JUT d max |/( e m p6»)|— £ (63) 
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Proof. Let 6 E S^ 1 . The map ip 9 : {z g C : \z\ < t} -> (0) defined 
by tpg (z) = zO is clearly holomorphic. Hence fg defined by fg (z) = f (z9) = 
f o (pg (z) is holomorphic. It follows that f k ,i (z) defined by 

fk.i{z)= f f{ze)Y k , l {6)d9 (64) 

is a holomorphic extension of fk.i to {z € C : \z\ < r}. Cauchy's inequality 
shows that for p = \z\ 

\fk,i(z)\ 2 < I \f{z0)?d6- f \Y ktl (8)\ 2 d6 

< oj d max \f (e lu p9)\ 2 . 

The Cauchy estimates \g^ (0)| < max| z | =p \g (z)\ for a holomorphic function 
<7 and the last estimate imply for < p < r 



dzm+k A.' (°) 



max |/(e m pg)| ■ ( fc+ J ! ( 65 ) 

Let us write f k ,i (z) = J2m=k ir^r/fcj (0) ■ z m for \z\ < t. It is known that 
r~ k ' fk.i (r) is an even function, hence we can write 

00 i l2m+k 

p kll {r 2 ) = r-* h , (r) = £ (fc + 2m) , dr2ro+fe /M (0) • 

m— 

Then for t = r 2 we have 

■^Pk.i (t) = J! (fc + 2 m)! (m-s)!dr 2 ™+ fe/M (0) ' ^ ' 
Now (1651) implies 



d s 

Ws Pk,l (t) 



1 °° 

< max 1/ (e m p9) I — - V 



tte[o,27rl,6>eS< 1 - 1 1 pfe+2s (m — s)! \P 2 

m! 4-rn—s " x y-m j 

(m-s)! 1 "~ dt s Z^m=0 1 "~ tit 8 1-t 



Since for \t\ < 1 we have V°° "l> i™- s = *1 V°° _ t m = d ° 

s! (1 — t) & 1 a straightforward computation yields (|63[) . ■ 

Now combining Theorem 1261 and the last Lemma, we obtain the final esti- 
mate. 

Theorem 29 Let < R < oo and Zei p be a pseudo-positive signed measure 
with support in Br satisfying the summability condition \5%fy and let be 
the polyharmonic Gaufl-Jacobi measure of order s. Then the following error 
estimate 



r— 2 00 a fc -i i-R 2 

E(f)< max 1/ HI EE -u / \Qi 

(p 2 — R 2 ) wec,\w\< P t~Qi~i p Ja 



i(t)\ (*) 



Zio/ds /or aZ/ functions f : Br — > C which possess a holomorphic extension to 
the complex ball B^ for r > R and for any p with R < p < r. 
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6 Polyharmonic Gaufi-Jacobi cubature in the 
annulus 



Let us imagine a function / which is say holomorphic in the ball Br with some 
singularities in the smaller ball B°, and one needs to find the integral of the 
function on the annulus A PiR with an estimate for the error of approximation. 
This example is a motivation to consider polyharmonic Gaufi-Jacobi cubatures 
in the annulus A PiR which generalize the construction of Section 

By the results in Section 0] the polyharmonic Gaufi-Jacobi measure of 
order s has support in A PtR and it is exact on the space of all polynomials f such 
that A 2s f = 0. In the present section we will seek polyharmonic Gaufi-Jacobi 
cubatures which are exact on the larger space 

PH 2s (A p , R ) = {feC (A p , R ) n C 2s (A° PtR ) : A 2s f (x) = for all x G A° R } 

where A° R is the interior of A P . R . Exactness with respect to PH 2s (A p . R ) is 
related to the expectation that the integrals of functions with singularities in 
the inner open ball < p} will be better approximated. 

It turns out that the problem can be solved in a way very similar to Theorem 
1201 The proof is so far based on the "generalized Gaufi-Jacobi quadratures for 
Chebyshev systems" which have been developed mainly by A. Markov, |.'S2I 
Chapter 4]. We restrict our discussion to the case of compact annulus which is 
less technical. Let us now formulate the result precisely: 

Theorem 30 Let < p < R < oo. Let p be a pseudo-positive signed measure 
with support in A PtR such that 

r- k dp k ,i(r) < oo. (66) 

k=0 1=1 " p 

Let s be a natural number. Then there exists unique pseudo-positive measure 
T ( 2s ) w ith support in A p jf such that 

(i) The cardinality of the support of each component measure Tt j , k £ No, 
I = 1,2, dfc, is < 2s. 

(ii) J f (x) dT<M (x)=Jf (x) dp (x) for all f e PH 2s (A PjR ) . 

Let us compare the result with Theorem [201 In the latter case we obtained 
a measure with support in A p R which is exact for all polynomials P with 
A 2s P = 0. The support of the component measure crf s ; has at most s points. 
In contrast, the component measure rj^" of the solution t^ 2s ^ has a support of 
cardinality < 2s which is twice bigger. This is caused by the fact that t^ 2 ^ is 
exact on the larger subspace PH 2s (A PjR ) . 

For the proof we will first need the representation of a polyharmonic function 
in the annulus which is somewhat more sophisticated than that of a polynomial 
as given in formula l|29|l . Let us introduce the operators 



d 2 n - 1 d k(n + k-2) 
L ^ : = d ^ + — d -r ~ (6?) 
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which may be written as 

1 d r .J r 1 11 

(68) 



+fc-i dr 



dr 



see e.g. (10.18) in Let L?|sdenote the 2s-th iterate of These operators 
are the radial part of the polyharmonic operators A 2s . A set of the 4s linearly 
independent solutions of the equation 

L 2 { k } f(r) = forr>0, (69) 

can be explicitly constructed; e.g. for k > 4s a set of solutions is given by 

r 2j - k j = 0,l,...,s-l; r^'+ fe , j = 0, 1, 2s - 1. (70) 

For k < 4s one has to be careful with multiplicities, and we refer to f° r an 
explicit description. 

We have the following result (see e.g. Theorem 10.39 in |31|): 

Proposition 31 Let h e PH 2s (A Pi r) . Then the Laplace-Fourier series 

oo d k 
fc=0 1=1 

converges absolutely and uniformly on compact subsets of A° R . The Laplace- 
Fourier coefficients hki are solutions of fffffi . 



In the following wc will mimic the proof of Theorem 1201 The measures a 



O) 

k.l 



in the proof of Thcorcm|2nihad the feature that they were exact on the solutions 
l,r 2 , ...,r 4s ~ 2 . Proposition 1311 shows that we need now quadratures which are 
exact on the solutions of (|69J) . This motivates to recall the theory of A. Markov 
and M. Krcin on quadratures for Chebyshev systems. 

Definition 32 Let uq,...,un be continuous functions on [a, b] . We say that 
they form a Chebyshev system of order N + 1 on [a, b] if every non-trivial linear 
combination Ylj=olj u j (*) nas °^ m °st N zeros on [a,b] , i.e. the determinant 

det ( Uj (t,))5 =0 (72) 

is not zero on [a, b] . The system u , Itjv is T + on [a, b] if det (v,j (ti)) i -_ > 
holds for every choice of tj € [a, b] with t < t\ < ... < t^. 

Note that the definition of a T + -system depends on the order of the functions 
Uq, ...,u N . 

The following theorem is the generalization of the GauB-Jacobi quadratures 
for Chebyshev systems (see Theorem 1.1 and Remark 1.2 in Chapter 4, and 
Theorem 1.4 in Chapter 2 in (321)- 
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Theorem 33 Let N = 2s — 1 for a natural number s and let a be a non- 
negative measure on [a,b] with cardinality of the support > s. Let the continuous 
functions uq, ujy be a Chebyshev system on the interval [a, b] and assume that 
un+i := £1 is a continuous function on [a, b] . If uq, un + i is a T + -system on 
[a, b] then there exists a unique measure with support of cardinality s such 
that 

rb fb 

Uj (t)da(t)= Uj{t)da^{t) for j = 0,1,..., N. (73) 



"3 VI VI — I <*3 
J a 

The support of is contained in the open interval (a, b). 

The measure is called in [221 the "lower chief representation" which is 
also very natural to be called Gaufi-Jacobi-Markov quadrature, and we will use 
this name further. 

A second major result in the Krein-Markov theory is the extremal property 
of the truncated moment problem due to Chebyshev, Markov and Stieltjes. 

Theorem 34 With the notations and assumptions of Theorem \3!ft let <jW be the 
Gaufi-Jacobi-Markov quadrature of a. The measure attains the minimum 
in the problem 

n(t)dv(t) (74) 
where v ranges over all non-negative measures v such that 

b rb 

Uj{i)dv(t)= / Uj(t)da(t) for j = 0, 1, TV. 

J a 

We return now to our case of polyharmonic cubatures on annuli. At first we 
note 

Proposition 35 Any linear independent system R s k j,j = 1,2,. ..,4s, of solu- 
tions of \b is a Chebyshev system of order 2s on every interval [a, b] with 
a > 0. For k > 2s the system of solutions in is a T + — system on [a, b] with 
a > 0. 

The proof follows from the results in Section II. 5 and Theorem II. 5. 2 in (35] 
and uses the representation I|tj8fl ■ In view of the Krein-Markov theory we need 
the following stronger result: 

Proposition 36 Let k E N , I = 1,2, ...,a/- be fixed, with k > 4s, and define 
u is = Q = 1. If we denote by Uq, u^ s -i the system of solutions in ^7u\ ) then 
uq, U4 S is a T + — system on [a, b) for a > 0. 

Proof. By the example in |32l Chapter II, Section 2.1/c)] the system 
{e^- k \ j = 0,1,. ..,28-1; e°; e^+ k \ j = 0,1, ...,2s - l} 
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is a T + — system since the numbers 2j — k and 2j + k are all different due to 
A; > 4s. Then the reordered system 

[e^~ k \ j = 0,l,...,2 S -l; e^+ fe >, j = 0, 1, 2s - 1; e } 

has the same determinant sign of ((72(1 as the above. By a change of the variable 
r — e*, one concludes that the system uq, ...,it4 S is a T + — system. ■ 
Now we are prepared to make the proof. 

Proof of Theorem 1301 Fix a pair of indices (k,l) with k £ No, / = 
1,2, ...,afc, and let fik,i be the component measure. If the support of fik,l has 
less than 2s points, put t^ 2 ^ = (ik,i- Assume now it has more than 2s points. 
Let Uq, Ui s -i be the system H70I) . By Markov-Krein's Theorem 1331 applied 

(2s) 

to a := fj,k,i there exists a GauB-Jacobi-Markov measure ; with support in 
(p, i?), and its support has cardinality 2s. An essential point is to prove that (at 
least) for sufficiently large k one has 

[ R ldT^(r) < j ldfj,k,i (r) ■ (75) 

For k > 4s this follows immediately from Markov-Krein's Theorem l34l bv means 
of Proposition 1361 

Further we proceed as in the proof of Theorem [20| We want to define a 
functional T^ 2s ) on C {A p , R ) by putting 

OO O-k pR 

^(/):=EE fkAr)drff(r) (76) 
k=0 1=1 Jp 

for / 6 C (Ap^n) where fk,i are its Laplace-Fourier coefficients. Indeed, by using 
the standard estimate for the Laplace-Fourier coefficients the inequality 

oo a k „R 

T<M(f) <C max I/WIEE / ^ (r) 

pS|2;|S k=0l=l J P 

is easily established for all / € C (A Pi #). Now with (|75|l and our assumption 
1)66(1 it follows that T^ 2s ^ is well-defined. We may apply the Riesz representation 
theorem and obtain a representing measure, denoted by t^ 2s \ with support in 
A p r. Since the constant function is in C (^4 p ,_r) it is clear that t^ 2 ^ is a finite 
measure. Let us remark that due to our assumption 1(6*61) the identity 



/ oo a k „R 
f( X )d^ = j2J2 / /*. 
k=0l=l J P 



i (r) d/z fc ,i (r) 



holds for all / G C (^4 p ,,r) , since the right hand side defines a continuous func- 
tional on C (A p ^) which agrees with fi on the dense subspace C x (A Pj r). Due 
to the exactness property of all measures T^f 1 and the representation 1(76(1 . it 
follows that t( 2s ) satisfies ii) of Theorem EU for all / € Pi? 2s (A p , fl ). 
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The pseudo-positivity of r( 2s ) follows from Corollarvll 21 since T( 2s ) is clearly 
pseudo-positive definite. The uniqueness of t^ 2 ^ follows from the uniqueness of 
the Gaufi-Jacobi-Markov measure rj?^ as in the proof of Theorem [20| ■ 

7 Polyharmonic Gaufi— Jacobi cubature in the 
cylinder (periodic strip) 

The concept of pseudo-positivity which we have studied so far depends on the 
expansion of the polyharmonic functions in Laplace-Fourier series which uses 
the rotational symmetry of the ball and the annulus. The polyharmonic Gaufi- 
Jacobi cubature in the ball was defined respectively by the application of the 
GauB-Jacobi quadrature to the Laplace-Fourier coefficients. It is natural to 
extend this concept on expansions available in other domains with symmetries, 
and we will do so for the case of the cylinder (which may be considered also as 
a periodic strip), where the Fourier series is the natural expansion. 

Let now -co < a < b < oo. We consider functions / : [a, b] x R d_1 — > C 
depending oni= (t, y) in the strip [a, 6] x R d_1 which are 27r-periodic 10 with 
respect to the variable y, i.e. satisfy the equality 

f(t,y + 2vrfc) = / 0, y) for all i € R, y G M d ~\ k € Z d ~ l . 

Let us introduce the cylinder 

X = [a,b] x T^ 1 

where T d_1 = S d_1 is the d — 1 dimensional torus. One may interpret the 
space X also as a periodic strip. The space of 2~k— periodic in y functions 
coincides with the space C (X) . By PH S ([a, b] x T d_1 ) we denote the space 
of functions / : X — > C which are polyharmonic of order s on X, i.e. the 
space of functions / € X such that the corresponding 2ir— periodic in y function 
/ G C ([a, b] x K d_1 ) is polyharmonic of order s in [a, b] x R d_1 of the variables 
x. 

If we fix t £ [a, b] , the equality 

f(t,y)= fk(t)e i{k ' v) (77) 

is the Fourier series expansion of / € C (X) where the Fourier coefficients fk (t) , 
k S Z d_1 , of / (t, •) are given by 

(2?r) Jt^- 1 (2tt) Jo Jo 
(78) 

10 The case of non— periodic functions on the strip [a,b] X M d_1 is similar but needs more 
care. 
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In the following we want to construct polyharmonic GauB-Jacobi cuba- 
tures for measures u defined on R x T d_1 , or equivalently, denned on R d and 
27T— periodic with respect to the variable y. 

Next we introduce pseudo-positivity in this setting: 

Definition 37 A measure /i on R x T d_1 with support in the cylinder X is 
pseudo-positive (in order to avoid mixing with the notion of pseudo-positivity 
introduced in Section\^we will say sometimes pseudo-positive on the cylin- 
der X), if for each k £ Z d_1 and for each non-negative continuous function 
h : R — > [0, oo) with compact support the inequality 

[ h{t)e l{k ' v) d^{t,y) > 
Jx 

holds. 

For every k £ Z d_1 , one may apply the Riesz representation theorem to 
prove the existence of a unique non-negative measure fik on R such that 

h (t) dn k (t) = / h (t) e^ dfi (t,y) (79) 
Jx 

holds for any h £ C c (R) . If we assume that [i has support in X = [a, b] x 
T d_1 with -oo < a < b < oo then it is clear that J7SJ| holds for all h € C (R) . 

Now we are going to prove the existence of polyharmonic GauB-Jacobi cu- 
bature for the case of the cylinder. 

Theorem 38 Let — oo < a < b < oo. Let n be a finite, signed measure with 
support in the cylinder X = [a, b] x T d_1 which is pseudo-positive on X . Suppose 
that 




<i/i (t, y) < oo. 



Then for each natural number s there exists a unique finite signed measure a^ 2 '^ 
with support in X such that 

(i) The support of each component measure ,k £ 7L d ~ x is in [a, b] and 
has cardinality < 2s. 

(ii) J fdu = J fdo-W for all functions f £ PH 2s (X) . 
(Hi) a^ 2s ^ is pseudo-positive on the cylinder X . 

Proof. Let us give a characterization of the space PH 2s (X) . If / £ 
PH 2s (X) then / is representable in the Fourier series H77J1 where for every 
k £ Z™ _1 the function ft (t) is a C°°-solution to the equation 

(«) = „, 
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cf. Theorem 9.3 in All solutions of the latter equation are linear combina- 
tions of the following functions 

Uj (i) = ^e~ |A:| * for j = 0, 1, 2s - 1, 
u 2s+j (t) = Ve^ for j = 0, 1, ...,2s - 1. 

Further define the function Ui s = 1. Then the system of functions uq, ■■■,M4s is 
a T + system - this follows from example c) in Chapter 2.2 of |32) . 

Let /ik be the non- negative measure defined by l(?9"jl . If the support of fit has 

(2s) 

less or equal than 2s points we define a k :— /ij,. If the support of fi^ has more 
than 2s points, there exists according to Theorem 1331 a non-negative measure 
cr^ such that 



b 



Uj (t) dcr[ 2s) 0) = / Uj (t) dn k (t) for j = 0, 4s - 1, 

J a 

and the support of <J^ S ^ has < 2s points and lies in [a, b] . From the Krein- 
Markov Theorem 1341 it follows 

ldcr k (t) < { Id^u (t) . (80) 



Following the usual scheme, we want to define the functional T^ 2 ^ on 
C (X) by putting 



t {2s] (/) - — U E / 'mWP 



We have to show that the functional T^ 2 ^ is well-defined. Note that 

T^(f)<-^r E max|/*(*)l /l*Tfc(*). (81) 

Moreover it is clear that 

max|/fc(t)|< max max 

t£[a,b] t£[a,b] yET 1 - 1 

By (jHOJ we obtain 

T^Hf) < max max \f(t,y)\ T [ ld/i* (t) , 

Thus the functional T^ 2s ) (/) is well-defined on C (X) , and by the Riesz rep- 
resentation theorem there exists a signed representing measure with support in 
X, denoted by a^ 2s \ Arguments similar to those presented at the end of the 
proof of Theorem [3U| show that a 1 * 2 ^ is pseudo-positive, unique and that the 
exactness property (ii) is satisfied. The details are omitted. ■ 
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8 Examples and miscellaneous results 

In this section we provide some examples and results on pseudo-positive mea- 
sures which throw more light on these new notions. 

8.1 The univariate case 

It is instructive to consider the univariate case of our theory: then d = 1, 
S° = { — 1,1}, and the normalized measure is u>o(6) = \ for all 9 E S°. The 
harmonic polynomials are the linear functions, their basis are the two functions 
defined by Yq (x) — 1 and Y\ (x) = x for all x G M. The following is now 
immediate from the definitions: 

Proposition 39 Let d = 1. A functional T : C [x] — > C is pseudo-positive 
definite if and only if T [p* (x 2 ) p (x 2 )) > and T (xp* (a; 2 ) p (x 2 )) > for all 
p€C[x}. 

It follows from the last proposition that a Stieltjes moment sequence is always 
pseudo-positive definite; by definition the functional T : C [x] — > C has the 
stronger property that T (q* (x) q (x)) > and T (xq* (x) q(x)) > for all q £ 
C [x] . Below we give an example of a pseudo-positive definite functional which is 
not positive definite, in particular it does not define a Stieltjes moment sequence. 

As pointed out in |451 Chapter 4.1], the Laplace-Fourier expansion of / is 
given by 



are the usual even and odd functions. 

Example 40 Let a be a non-negative finite measure on the interval [a, b] with 
a > 0. Then the functional T : C [x] — > C defined by 



is pseudo-positive definite but not positive definite. 

Proof. Let fi be as above. Then T (/) = 2 J fi (r) da, so T is pseudo- 
positive definite by Proposition |SJ Since T(l) = and T ^ it is clear that T 
is not positive definite. ■ 






f(-x)da 
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8.2 A criterion for pseudo-positivity 

The following is a simple criterion for pseudo-positivity: 

Proposition 41 Let /i be a signed moment measure on M. d . Assume that /i has 
a density w (x) with respect to the Lebesgue measure dx such that 9 i — ► w (r&) 
is in L 2 (§ d_1 ) for each r > 0. If the Laplace- Fourier coefficients of w, 

w k ,i (r) := f w {r9) Y k<l (9) d9 

are non-negative then fj, is pseudo-positive and 

d m (r) = r fc+d - V,* (r) , (82) 

r- k dii k:l (r) = / w k ,i (r) ■ r d ~ l dr (83) 
o Jo 

if the last integral exists. The measures fik,l are defined by means of equality 

Proof. Since /i has a density w (x) we can use polar coordinates to obtain 
for / G C pol (R d ) 

fdfi= [ f{x)w(x)dx= [ [ f (r6)w(r6)r d - l d6dr. (84) 
JWL d Jo Js d - 1 

For any h G C po i [0, oo) we put / (x) — h (\x\) Y^.i (x) , then we obtain 

h(\x\)Y Kl (x)d l i= f f h{r)r k+d - 1 Y Kl {e)w{re)d9dr. (85) 
Jo jV- 1 

Since 9\ — > w (r9) is in L 2 (S^ 1 ) , we know that w kJ (r) = J gd _ 1 w (rO) Y kJ (9) d6. 
Hence, by the definition of fJ>k,l, we obtain 

I* /"OC 

h(r)dnk,i~ h(\x\)Y k j(x)dn= h(r)w k , l (r)r k+d - 1 dr. (86) 

JTSL d JO 

Thus the measure fi is pseudo-positive, and l|82|l follows. Let us prove l|83() : we 
define the cut-off functions h m € C po i [0, oo) such that h m (t) = t~ k for t > 1/m 
and such that h m < h m +i - Now use H86(l and the monotone convergence theorem 
to obtain (jHSJ- ■ 

The next example addresses the question of whether there is a relationship 
among the supports of the GauB-Jacobi quadratures a k} i in Theorem 1201 

Proposition 42 Let s be a natural number. Then there exists a pseudo-positive 

(s) 

measure /i with support in the unit ball such that the component measures a k [ of 
the polyharmonic Gaufi-Jacobi cubature in Theorem \2(A of order s have identical 
supports. 
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Proof. Let < a < b and consider the density 



oo / c _|_ c /_x -j 

fc=i 

where l[ a ,6] is the indicator function of [a, b] . Then w induces a pseudo-positive 

measure and J pdfik.i = — p — J " 2 p(t)dt according to l|86[l . It follows that 
for all k £ N the orthogonal polynomials of degree s associated with /j,k,i are 
identical up to a factor. Hence the supports of the measures <J^\ are identical 
for aU k £ N. ■ 



8.3 The two— dimensional case 

Let us consider the case d = 2, and take the usual orthonormal basis of solid 
harmonics, defined by Yq (e 1 *) = and 

y M (re lt ) = -j=r k cos kt and Y K2 (re lt ) = -^=r k sin kt for k £ N. (87) 

We define a density to( a ) : R" — > [0, oo), depending on parameter a > 0, by 

j (q) (re lt ) := (1 - r Q ) P (re lt ) for < r < 1 
. ; («) ( re «) _ o for r > 1; 



here the function P (re lt ) is the Poisson kernel for < r < 1 given by (see e.g. 
5.1.16 in 2, p. 243]) 

1 - r 2 °° 
P (re lt ) := T = 1 + V 2r fc cos kt. (88) 



fc=l 



By Proposition I411 the measure d\x a := (x) dx is pseudo-positive. For 
k > 0, by JS3JI and (JH3 we obtain 



.fc+i 



r 



(1 - r a )dr = 



2^c 



(A + 2) {a + k + 2) ' 



It follows that w^ a > (x) dx satisfies the summability condition 111(1 . so we can 
apply our cubature formula to this kind of measures. 

On the other hand, there exist pseudo-positive measures which do not satisfy 
the summability condition Ijllfl: 

Proposition 43 Let w (re lt ) :— P (re lt ) for < r < 1 and w (re lt ) := for 
r > 1 where P (x) is given by ]8fy) . Then d^i := u> (x) <ix is a pseudo-positive, 
non-negative moment measure which does not satisfy the summability condition 
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Proof. It follows from for k > 1 

J r- k dfi k s = J™ w kA (r) ■ r d ~ l dr = 2yfr £ r k+1 dr = 

so we see that the summability condition is not fulfilled. ■ 

Next we compute explicitly the error in Section [S] for the function (x) 
with a — 2. 

Theorem 44 Let dfi :— w^ 2 ' (x) dx, and let crM be the polyharmonic Gaufi- 
Jacobi measure of order s. Then for every f € C 2s (R d ) the error E s (f) 
| f (x) dfj, — f f (x) da 1 -^ can be estimated by 



(2a)! 



oo _j2s 
k=0 v<ti< 



s!(s + fc + l)!(s + l)!(5 + fc)! 
(2s + fc + l)! (2s + fc + 2)! 



Proof. For k > we obtain by Ij86(l and (|88|l the equality 
p (t) dfi kA = 2^ / P (r) (1 - r 2 ) r 2k+1 dr, 



and clearly = 0. Let ip be the usual one defined by ip (t) = t 2 . Then fif 1 
can be computed by 



p (t) d^ x = / p (i 2 ) dn ktl = 2^ p (r 2 ) (1 - r 2 ) r 2k+1 dr. (89) 



According to Theorem 1261 we have only to compute J Q Q s kl (t) dfi k[ where 
Q s k j (t) is the s-th orthogonal polynomial with leading coefficient 1. The substi- 
tution t = r 2 in the integral JS^J yields J,, 1 p (t) dfxf t = J 1 p (t) (1 — f ) t k dt. 
The substitution < := | (x + 1) and dt — |dx shows that 



f P W <i = ^1 / / Q ( x + !)) (1 " • (1 + *)* 



dx. 



Let Pg"'^ (x) be the Jacobi polynomial of degree s (see p. 30]) normal- 
ized with p^ a ''^ (l) — ( S ^ Q ). They are orthogonal with respect to the measure 
dvja^p := (1 — x) a (1 + x) dx for a > —1, (3 > —1. It is known that the leading 
coefficient k s of P$ a>f3) (t) is equal to k s = 2- s ( 2s +^+^). Further 



dw 



a,0 



2 a +f 3+1 r(s + a + l)T(s + /3 + l) 
2s + a + (3 + 1 s\T(s + a + (3+ 1) ' 



Define P s (Q:/3) (i) := P s (Q:/3) (2t - 1) for < t < 1. Then P^ k) are orthogonal 
polynomials for 1 since 

/ P^P^dd, 



2 k+2 



{x) p£,*) (i _ X ) (i + j)* (go) 
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The leading coefficient fc s of pj 1 '^ (t) = pj 1 '^ (2t— 1) is equal to 2 s k Sl so 

k. 



k s = ( 2s+ s fe+1 ). Thus Q s kl (t) := -i-P s (1 ' fc) and we obtain from P7|l 



1 V^t 

£ 22fc+2 - < 



P^ hk) (x) (1 - x) (1 + x) k dx. 



Now this is equal to 

'2s + k + 1\ _1 2 fe+1+1 r (s + 2) r (s + k + 1) 



giving 



for large fc. 



2 fe + 2 2s + fc + 2 sir (s + fc + 2) 



2 „ _ s!(s + fc + l)!(s + l)!(s + fc)! 



1 



(2s + fc + l)! (2s + fc + 2)! fc 2s + 2 



8.4 The summability condition 

We are now turning back to the general situation. The next result shows that 
the spectrum of the component measures a k ,i is contained in the spectrum of 
the representation measure \x. 

Theorem 45 Leta k i be non-negative measures on [0,oo). If the functional T : 
<C[xi, X2, Xd] — ► C defined by 113 3p possesses a representing moment measure 
H with compact support then 



0~k 



,*({m 2 })< max in,, wi-i^-im^N 2 ^- 1 



for any x E M. d where \fi\ is the total variation and \x\ 2 E> d 1 = {|a;| 2 # : 9 £ 

the support of /i be contained in Br. Let xq € M. d be given. For 
" ■ i-.i with p ^|xq| 2 ^ =1 we have 



Proof. Let 

every univariate polynomial p 



,i({\xo\ 2 }) < r p{r 2 )da kJ < f p(\x\ 2 )Y k>l (x) 
v 'Jo Jm d 

< max |Y fc ,,(0)| / P (\x\ 2 ) \x\ k d\fi\. 



2 J = 1 which converges 

2 



Now choose a sequence of polynomials p m with p m (\xq\ 

on [0, i?] to the function / defined by / (j£o| 2 ) =1 an d / (t) = for t ^= \xq 
Since |/z| has support in Br Lebesgue's convergence theorem shows that 

\Yk,i((>)\ I \f{x)\\x\ k d\n\. 



<7k, l 



({kol 2 } 



< max 
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The last implies our statement. ■ 

The following result shows that the summability condition is sometimes 
equivalent to the existence of a pseudo-positive representing measure: 

Corollary 46 Let d = 2. Let ak,i be non-negative measures on [0, oo) and as- 
sume that they have disjoint and at most countable supports. If the functional 
T : C[xi,X2] — > C defined by possesses a representing moment measure 
with compact support then 



00 a k p 00 

VV / r- k da kd (r) < 00. 

r._n JO 



k=0 1=1 



Proof. Let be the support set of a k ,i- The last theorem shows that 
o~k,i ({0}) = 0, hence ^ Moreover it tells us that 

/>oo 

/ r- k da Kl (r) < max \Y k>l (6)\ ■ V | M | (rS^ 1 ) . 
Jo ees^ 1 7^ 

Since d = 2 we know that max# g gd-i \Yk,i (0)\ < 1. Hence 

00 a k „oo 00 a k 

Y: E / r- h *w «<EEE ^ Ia»i (r") 

fe=0 i=i ■ J ° k=0 1=1 res fci! 

where the last inequality follows from the fact that E^,; are pairwise disjoint. ■ 
Recall that the converse of the last theorem holds under the additional as- 
sumption that the supports of all Uh,\ are contained in some interval [0, R] . 

Theorem 47 There exists afunctional T : C [x\, X2, xj] — * C which is 
pseudo-positive definite but does not possess a pseudo-positive representing mea- 
sure. 

Proof. Let a be a non-negative measure over [0, R] . Let / G C [sci, .., stj] 
and let be the Laplace-Fourier coefficients of /. By Proposition it is clear 
that 

T{f) := / /y(r)r^(r) 
Jo 

is pseudo-positive definite. We take now for tr the Dirac functional at r = 0. 
Suppose that T has a signed representing measure /x which is pseudo-positive. 
Then the component measure is non-negative, and it is defined by the 
equation 



h(r)diiu(r) := / h(\x\)Y 11 (x)d f i 

for any continuous function h : [0, 00) — > C with compact support. Take now 
h (r) = r 2 . Then by Proposition |H1 

f r 2 d^ n (r)= f \x\ 2 Y n (x) dp = T (|x| 2 Y n (*)) = 0. 
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It follows that /in has support {0} . On the other hand, if we take a sequence 
of functions h m € C c ([0, oo)) such that h m — > l{o}: then we obtain 

A*n ({°}) = lim / h m (\x\)Yu(x)dfi. 

But h m (\x\) Yu (x) converges to the zero-function, and Lebesgue's theorem 
shows that /in ({0}) = 0, so fin — 0. This is a contradiction since 

roc p pR 

/ ldnu(r)= Yu(x)dn = T{Y n )= lda(r) = l. 

JO JR™ JO 

The proof is complete. ■ 



9 Concluding Remarks 

One important feature of the polyharmonic Gaufi-Jacobi cubature which de- 
serves to be discussed is its numerical significance. At first glance, one may 
object that the cubature needs the knowledge of the Laplace-Fourier coeffi- 
cients of the function / : R" — > C, and these are based on integrals as well. 
However, if one works with polynomials, the Gaufi decomposition can be 
constructed by an efficient differentiation algorithm, see 0]. Decomposing the 
harmonic polynomials hj according to our fixed orthonormal basis Yk,i (x) one 
obtains from the expansion 

deg / a k 

f(x)=J2J^ Pk ,i(M 2 )Y k ,t(x), 

k=0 1=1 

where p kt i are uniquely determined univariate polynomials. For the Gaufi-Jacobi 
quadrature there exists suitable software to compute the weights , aj'' 
and the nodes r[ k ' l \ ri k ' 1 ^ of the measures cr^\. Then 

dog/ a k s 

T is) (f) :=£££a?Vi Of ) 2 )- 

k=0 1=1 3=1 

In practice one also has to bound the number of spherical harmonics Y k> i (x) 
in the formula. Our main results Theorem I2UI and Theorem 1221 show that by 
increasing the number k and the number s the algorithm remains stable. 

What concerns the polyharmonic Gaufi-Jacobi cubature in the annulus and 
the strip, we note that recently efficient algorithms for finding Gaufi-Jacobi- 
Markov quadratures for Chebyshev systems have been studied in 

A second point to be made clear, is the motivation why we choose the space 
of polyharmonic functions of order s as the exactness space for the multivariate 
generalization. One main reason is the fact that recently polyharmonic functions 
have shown to be an efficient tool in approximation theory and more generally, in 
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mathematical analysis, see e.g. [3], (2E1) |30| . |31| . and Another motivation 
stems from potential theory: two non-negative measures fi and v with compact 
supports are gravitationally equivalent if 



for all harmonic functions h defined on a neighborhood of the supports of the 
measures. If /i and v are gravitationally equivalent then they produce the same 
potential outside of their support. Graviequivalent measures are very important 
in inverse problems in Geophysics and Geodesy, and a new mathematical area 
has grown extensively during the last two decades or so under the title Quadra- 
ture Domains, see the comprehensive survey and references in |23| . as well as 



In analogy to (|91|l one could define two (generally speaking, signed) mea- 
sures fx and v as polyharmonically equivalent of order s if (|91Jl holds for all 
polyharmonic functions of order s in a neighborhood of their support. Similar 
notions of equivalence have been developed by L. Ehrenpreis m the form 
of a generalized balayage. To our knowledge the equivalence of two measures 
(and more generally, distributions) with respect to the solutions of an elliptic 
operator has been for the first time rigorously formulated and studied in the 
case of non-negativity in 021 ■ 

Let us remark that the polyharmonic GauB-Jacobi measure which we 
have introduced in the present paper is related to the concept of "mother body" 
(a non-negative measure v satisfying l|91|l for a given [i, and having minimal 
support) in the theory of Quadrature Domains, cf. [221, and Remark|23 as well 
as (23 and 
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